Stability of trapped Bose-Einstein condensates * 
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In three-dimensional trapped Bose-Einstein condensate (BEC), described by the time-dependent 
Gross-Pitaevskii-Ginzburg equation, we study the effect of initial conditions on stability using a 
Gaussian variational approach and exact numerical simulations. We also discuss the validity of the 
criterion for stability suggested by Vakhitov and Kolokolov. The maximum initial chirp (initial 
focusing defocusing of cloud) that can lead a stable condensate to collapse even before the number 
of atoms reaches its critical limit is obtained for several specific cases. When we consider two- and 
three-body nonlinear terms, with negative cubic and positive quintic terms, we have the conditions 
for the existence of two phases in the condensate. In this case, the magnitude of the oscillations 
between the two phases are studied considering sufficient large initial chirps. The occurrence of 
collapse in a BEC with repulsive two-body interaction is also shown to be possible. 
PACS numbers: 03.75.Fi, 32.80. Pj, 42.50.Md, 03.75.-b 



I. INTRODUCTION 

In trapped Bose-Einstein condensates(BEC), it is well- 
known the occurrence of collapse when the two-body in- 
teraction is attractive and if the number of atoms N ex- 
ceeds a critical value N,,, as in the case of atomic conden- 
sates with ^Li In this case, experiments with attrac- 
tive two-body interaction have been performed with 
results consistent with the limitation in the number of 
atoms and with the grow and collapse scenario. As the 
non-linear terms presented in the non-linear Schrodinger 
equation (NLSE) (cubic, quintic, and so on) are due to 
the first terms of an expansion of the effective many-body 
interaction in the mean-field approximation a study 
of stability of the equation that describes BEC is neces- 
sary, considering the observed few-body interactions that 
occur in the atomic gas. When the overall sign of the 
effective many-body interaction is negative, we have the 
conditions for the occurrence of collapse, with the stabil- 
ity of the condensate being restricted by the number of 
atoms. 

We should stress that the criteria for stability usually 
involves one parameter: the number of atoms. Here, in 
order to introduce a criterion for stability of a trapped 
condensate with two and three-body effective interac- 
tions, we introduce another parameter: a chirp, that is 
related to the initial focusing (defocusing) of the cloud. 
Without trap, a chirp parameter was previously intro- 
duced in Ref. [Q to study the stabihty of the NLSE; and 
with trap, the stability has been studied in Refs. [H^. 
One should also note that an oscillating condensate al- 
ways has a chirped wave function, where the chirp is pro- 
portional to the time variation of the width, which is a 
periodic function of time. Thus, a real condensate, for ex- 
ample in optical traps, always has nonzero initial chirp. 



that should be taken into account. 

Collapses for two and three dimension {2D and iD) 
geometries of the NLSE with trap potential has been in- 
vestigated in Refs. [|l]j8), where the corresponding critical 
numbers were obtained |^,|o| . These works are extensions 
of the Zakharov and co-workers theory of collapse [pd|-p[ 
to the case when an additional trap potential is included. 
The trap potential helps to prevent the collapsing pro- 
cess, which cannot be avoided if the number of atoms 
exceeds the critical limit. In Refs. [||,^, considering the 
moments method in the analysis of the equation for the 
average squared value of the width, (a^), it was empha- 
sized the influence of the corresponding initial condition 
on the collapse. 

Variational approaches have also been used by many 
authors for the analysis of stability of NLSE, and proved 
to be useful in many aspects (see Ref. Q, for example). 
The variational approach (VA), using a simple Gaussian 
ansatz, gives a reasonable description of the conditions 
for the collapse and also an approximate value for the 
critical number of atoms. It is interesting to apply the 
time dependent variational approach to describe the dy- 
namics of BEC with attractive two-body interactions and 
in particular to study the influence of the initial condi- 
tions (using a chirp constant, b, in the wave-function) on 
the stability of BEC. This can be important for the con- 
densation in optical traps and in strongly inhomogeneous 
traps. Also recently, in Ref. , the stability analysis was 
performed, using both VA and exact calculations, consid- 
ering that the condensate can unstabilize by tunneling, 
due to quantum fluctuations. 

When the two-body interaction term is switched off, for 
example by tunning an external magnetic field |p^ , the 
three-body interaction term can play an important role. 
Thus, it is interesting to study the collapse conditions in 
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the model with trap potential and an attractive three- 
body interaction (quintic) term. Recently, in Ref. [Q, 
it was investigated the collapse in one dimensional {ID) 
model with attractive three-body interaction. The trap 
was narrow and induced the radiative loss of atoms, re- 
ducing the number below the critical value, which made 
an arrest of the collapse possible. The two-dimensional 
{2D) case, with cubic nonlinearity (two-body interac- 
tions), may be useful in the modulation theory around 
the Townes soliton Then it is possible to estimate 

the role played by the continuum. 

In the present work, we first consider a review of 
the stability criteria, without application of initial chirp. 
Next, we analyze the role of an initial chirp in the wave- 
function of stable condensate with attractive two-body 
nonlinear term, which can be useful for actual experi- 
mental analysis in BEC. A generalization of this study 
of stability is considered with the inclusion of three-body 
interaction term in the NLSE. Such study can be relevant 
in the perspective of atomic systems with enhanced three- 
body effects, that can happen as the two-body scattering 
length is altered p7| . 

When a positive three-body interaction term is in- 
cluded in a trapped three-dimensional (3-D) NLSE with 
negative two-body interaction, it is already known the 
possibility of extending the region of stability of the equa- 
tion to larger number of atoms, with the occurrence of 
two stable phases in the condensate 21 1 . In Ref. |^ , it 
was calculated the frequency of the collective excitations 
in denser (liquid) and dilute (gas) phases of the conden- 
sate. Here, we also investigate a bifurcation phenomena 
related with the time variation of the width parameter of 
the theory, by switching the oscillations from one phase 
to the other. The present study is done with the aid of 
a time-dependent variational approach for the collapse of 
the atomic cloud in 31? model with harmonic trap po- 
tential. We compare the predictions of the VA with ex- 
act numerical simulations of partial differential equation 
(PDE). As detailed in our conclusions, the Gaussian VA 
gives a good estimate for the observables in the region 
of stability, and it starts to deviate from the exact re- 
sults when the system is close to the collapse conditions 
and particularly for the unstable solutions (maxima) of 
the total energy. For this regions, where the Gaussian 
VA fails, one should improve the ansatz or take the re- 
sults as a qualitative picture to guide the exact numerical 
calculations. 

In cases without traps, the NLSE was previously stud- 
ied by Vakhitov and Kolokolov (VK) where a cri- 
terion for stability was settled. This and other crite- 
ria for stability in nonlinear systems have been studied 
and extensively used by many authors (see (2^Jl^ and 
references therein). The VK criterion was recently de- 
tailed in Ref. When a trap potential was added to 
a NLSE with a negative cubic term, an alternative sta- 
bility criterion was derived in Ref. 0, after considering 
similarities between BEC atomic systems and compact 
object like neutron stars. Recently, the VK criterion for 
trapped BEC with cubic term was formally demonstrated 



in Ref. Q. 

In recent numerical studies one can also ob- 

serve that the VK criterion is not generally valid; it can- 
not be extended, for example, to the case of a trapped 
system with both attractive cubic and repulsive quintic 
terms. The solutions can be stable irrespectively to the 
sign of the derivative of the eigenvalues of the NLSE with 
respect to the number of atoms, dfi/dN. When 173 is pos- 
itive, and for large enough number os atoms, the NLSE 
is stable and dfi/dN can be positive (see, for example, 
Ref. Q), implying that the criterion cannot be extended 
to trapped NLSE with two and three-body terms. Here, 
the nonvalidity of the VK criterion is clearly verified in 
case we have a positive harmonic trap, with two and 
three-body nonlinear terms with opposite signs. 

We extend the study of the stability of the NLSE 
with cubic (two-body) and quintic (three-body) terms, 
through the analysis of chirp response, and through the 
frequency of collective excitations. By applying the per- 
turbation techniques to the theory of nonlinear oscilla- 
tions psj l , we also derive the frequency of weakly nonlin- 
ear collective excitations. This study extend numerical 
calculations performed in Ref. for the frequencies of 
collective excitations. 

One should also note that such study is of actual in- 
terest, as recently in Ref. |26| is was reported obser- 
vations of nonlinear oscillations in BEC of a gas with 
rubidium atoms. The present study can be of interest 
not only to atomic BEC but also to other branches of 
physics and mathematics, where non-linear effective in- 
teractions are added to a trapped potential in the time 
dependent Schrodinger equation, such as in optics and 
soliton physics. 

In Section II, it is given a description of the model, 
using a variational approach for a system with D dimen- 
sions. In Section III, we derive an expression for the chirp 
parameter, that is considered in our analysis of stability. 
The maximum initial chirp to keep the system oscillating 
in the same phase, around a minimum of the energy, with- 
out collapsing or without phase-transition, is obtained in 
this section by considering both variational and exact nu- 
merical calculations. In Section IV, we follow the study 
of Section HI, and also obtain the frequencies for the lin- 
ear and nonlinear oscillations. In this section we consider 
the presence of both, cubic and quintic terms, in several 
different configurations. The VK criterion of stability is 
discussed in section V. Finally, in Section VI, we present 
our concluding remarks. 



II. DESCRIPTION OF THE MODEL 

In this section we will consider the dynamics described 
by the NLSE with harmonic potential and with cubic and 
quintic terms. Using a variational approach for a system 
with D dimensions, we derive from the Lagrangian the 
expression for the anharmonic potential considering the 
mechanical analogy. 

The time-dependent Gross-Pitaevskii-Ginzburg (GPG) 
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equation that will be considered in the present approach, 
is given by 



2ii/;t 



-Ai/; + l^^ + A2|V^|V + A3|VlV- 



(1) 



In this equation and in the following, the explicit space 
and time dependences of the variables and parameters 
are implicit, unless it is necessary or convenient for clar- 
ity. The time derivatives will be denoted by indices t. 
In Eq. (0) we are assuming dimensionless variables: the 
unit of energy is fiuo/2\ the unit of length is y/Ti/ (mw); 
and the unit of time is \/uj. V = V{r) is a static trap 
potential, that we assume the harmonic oscillator with 
spherical symmetry and given by = 7'^, in the present 
units. A2 and A3 are the parameters of the two- and 
three-body interactions, which in general can be complex 
quantities. The imaginary parts of A2 and A3 describe, 
respectively, the effects of inelastic two and three-body 
collisions on the dynamics of BEC. In the present paper, 
we are not considering dissipative terms, and such cubic 
and quintic parameters are real. In order to compare with 
the formalism given in Ref. A2 is proportional to the 
two-body scattering length a^c and given by A2 = SttOsc. 
The chemical potential /i is given by the eigenvalue so- 



lutions of Eq. (|l|), with ijj{r, t) 



if{r): 



flip = -Aip + r^ip + X2\ip\'^ip + X^lifl'^if. 

The Lagrangian density corresponding to Eq. 
given by 



(2) 



IS 



C = 2 Imirt^j) - |V^|2 - r^l^p - _ ^|^|6. (3) 



2V2^ 3V3^ 



where 



6 J Ldt = 0, L = j £{r,t)dr. 



, (6) 



(7) 



Analyzing the corresponding Euler-Lagrange equations, 
dL d dL 



drji dt dfji 



0, 



(8) 



where rji refer to the variational parameters {A{t)^ o,{t), 
b{t) and 4>{t)), we obtain 



and 



at = ah] (9) 



A27Va-(^+2) 2A3iV2a-(2^+2) 



2(27r)(^/2) 



(10) 



Equation (j^ expresses the chirp parameter b via initial 
focusing (defocusing), (aj), of the wave- function. Com- 
bining the Eqs. and ([lO|), we have 



P 



att ^ ^ - a 



j(D+l) „(2Z5+1)' 



(11) 



A. Gaussian variational approach 

To analyze the dynamics of the BEC under two and 
three-body interactions, it is convenient to follow the 
variational approach developed in Refs. This ap- 

proach was successfully employed, recently, in for 
BEC with two-body interaction. We choose the simple 
Gaussian ansatz 



t) = A{t) exp 



.b{t)r^ 



(4) 



where A{t) is the amplitude, a{t) is the width, (f){t) is the 
linear phase of the condensate, and b{t) is the "chirp" 
parameter previously discussed. 

Without dissipative terms, the normalization of the 
wave-function is conserved and given by the number of 
particles N . The mean-square-radius and the normaliza- 
tion, for a system with D dimensions, are 



^ 2 , 
— a and 
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N = A^{^/^aY 



constant. 



(5) 



For a system with D dimensions and radial symmetry, 
the averaged Lagrangian expression is given by 



P 



Q 



X2N 

2^(2^^ 



(2)(^+3) A3 ^2 



inV3f (3)(^/^+i)A 



(12) 



Comparing the result of au with _D = 2, from Eq. (pT|), 
with Ref. psj , we observe that the resulting equation 
obtained by the variational approach coincides with the 
one obtained by the moments method |^9[] . The equa- 
tion is also close to the one obtained by the modula- 
tion method for the 2D Townes soliton ||l^. In case of 
D = 3, expressing the above definitions by the parame- 
ters given in Ref. ||2ll| , where A2 = SttOsc in our dimen- 
sionless units (flsc is the two-body scattering length), and 
N = n/{2V2\asc\): 



and Q = ^33. 



2V^ 



(13) 



The anharmonic potential is derived using the mechanical 
analogy. By using explicitly the dimensions in (||), 



R 




(14) 
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a result that can be identified with 

dUp{R) fnuo\ lY rmuj\^/2 dU 



OR 




da 



(16) 



In such extreme positions, we obtain 



Dai 



2a? 



{2-D) 



Q 



2 a 



2D 



{l~D), (23) 



where 



Up{R)^U{a) — . 



(17) 



The total energy for the BEC system with N particles is 
given by 



Et = NUp{R) - NU{a)^ = E{a)^. 

So, for the dimensionless quantities we obtain 
dU 



da 



= -{D)att. 



(18) 



(19) 



With (|1 Ij) , we have the anharmonic potential, in dimen- 
sionless units: 



rr, ^ D f ^ 1\ P Q 



(20) 



Correspondingly, from Eq. (y), the chemical potential is 
given by 



, ^ D f l\ 2P 3Q 



(21) 



The asymptotic limit a — s- oo, for C/(a), is the same 
of the oscillator [[/(oo) — s- oo], and is not determined by 
the parameters P and Q. So, the only limiting condition 
that is strongly affected by these parameters is the limit 
a = 0. When Q ^ 0, the three-body term is dominant 
in this limit: C/(0) w Q/{2a^^). Now, two cases have to 
be considered: Q < give us U{0) — > — oo; and Q > 
will give us U{0) +oo. As the other end is fixed by 
the oscillator condition [U{oo) — > -l-oo], in the first case 
(Q < 0) we cannot obtain more than one maximum and 
one minimum for finite a, that will depend on the sign 
and relative value of the parameter P. The case Q > 
can be very interesting, if P < 0: in both ends U{a) goes 
to -|-oo, such that in between we can have two minima 
and one maximum for the total energy. The two minima 
will represent two possible phases, and a rich dynamics 
can be described by the NLSE. From this analysis, we 
should note that phase transition is possible only in case 
Q > with P < 0. When Q < 0, there is only a single 
phase. 



B. Critical parameters for stability 



^{a,)^Dal + ^{i-D) + ^{3-D). (24) 

In Table I we resume the particular cases, from the 
above three equations, where exact analytical solutions 
can be found for the extremes of U and the correspond- 
ing critical parameters for stability. We note that, once 
a solution for a system with D dimensions and P ~ is 
known, the same solution can be used for the roots of a 
system with 2D dimensions and Q = 0, by just exchang- 
ing P and Q. When D = 2, the result Pc = —1 for the 
particular case when the quintic term is zero (Q — 0) is 
well known, and corresponds to = 2 in Rcf. ||30|| . The 
solutions will be unstable for P — — \P\ with \P\ > 1. 

Actually, one of the interesting cases occurs when the 
cubic term is zero (P = 0), for D = 1 and D = 3, in 
view of the possibility of altering the two-body scattering 
length by an external magnetic field . If this condition 
is realized (P = 0), the effect of the three-body non- linear 
term in the mean-field approximation of the condensate 
will be enhanced. 

TABLE I. For lower dimensions, we show the results of 
particular cases where exact analytical solutions can be found. 
as are the real and positive roots of the mean-square-radius 
that extremize the total energy [Et = NU{huj/2)]; Pc and Qc 
are the critical parameters for stability of the NLSE in just 
one local minimum. For D — 3, when P and Q are non-zero, 
we show just one critical set of parameters for the existence 
of two phases, f/c and fic are the corresponding critical values 
of U{a) and the chemical potential jJ.{a). 
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The relevant extremes of U{a) are given by the real 
and positive roots of 



In case such three-body (quintic) term is negative, as 
shown in Table I, the present variational ansatz gives 
analytical estimates for the critical parameters, and for 
physical quantities as the mean-square-radius, energy and 
chemical potentials. In Section IV we study in more de- 
tail this case. 

Other particular cases occur when phase transitions are 
possible. As explained before, these situations can only 
happen if both parameters P and Q are non-zero and 
Q > with P < 0. The absolute collapse is not pos- 
sible in the condensate; but, one can obtain transitions 
between two-phases. This implies that, for a system with 
two-body attractive interaction and three-body repulsive 
one, as we increase the number of atoms in the conden- 
sate (|P| increases), we can reach a critical limit where 
only one phase (stable) remains. In other words, for a 
particular value oi P ^ Pc, the position a of one of the 
minima of U coincides with the maximum, and we have 
an inflection point of U . 

As an example, in Table I we show one particular 
case where we have phase transition, and also analyt- 
ical solution for one of the roots of Eq. (|2^), when 
D — i. We can observe that this condition is realized 
if > P > = -(l/8)(i/'*). When Q = P^ the criti- 
cal limit occurs for Pc = —(1/8)'^^/'*-' and we have three- 
roots for as, given by a^i = 0^2 = \Pc\ = 0.5946 and 
asz =0.4696. 
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FIG. 1. For a stable BEG with A2 < 0, A3 > and 
D = 3, it is shown a variational plot of the potential en- 
ergy U{a) (upper frame), in units of {Tiuj)/2, as a function of 
the width a; and the corresponding phase plane defined by 
{at, a). The parameters are dimensionless, with P = —0.5495 
and Q =0.0099. a is related to the mean-square-radius as 
(r^) = 3aV2. 

A typical plot of the potential U{a) as a function of the 



width a and the corresponding phase portrait defined in 
the plane (at, a) for negative A2 and positive A3 is given 
in Fig. 1 for the case D = 3 of a stable BEG. The param- 
eters for the attractive cubic and repulsive quintic terms 
used in the figure are such that n =1.948 and =0.016, 
corresponding to P = —0.5495 and Q =0.0099. The num- 
ber of atoms corresponds to the situation in which we 
have two minima with the same value for U{a). From 
the variational expressions given in this section and our 
mechanical analogy, one should observe that, in the pic- 
torial example of Fig. 1, if we have a solution located at 
the right minimum (for example) it cannot migrate to the 
left minimum, unless the chirp parameter corresponds to 
an energy greater than the difference between the right 
minimum and the maximum that is in between. 



III. ANALYSIS OF THE OSCILLATIONS OF THE 
BEC IN THREE DIMENSIONS 

The different dynamical regimes in the condensate os- 
cillations can be described by the cross sections of the 
curve of U{a) with the levels of the effective total energy. 
Following our mechanical analogy, we have 



H{Pr,R) = ii + [/p(p) 
Zm 



(25) 



where Pf? is the momentum conjugate to R, given in (p^. 
In our dimensionless variables, the corresponding initial 
Hamiltonian is given by 



Ho 



U{ao) 



C/(ao), (26) 



where Eq. (|9|) was used. Thus, by varying the initial con- 
ditions ao , bo we alter Ho ■ 

This case was recently analyzed using the moments 
methods in Refs. ||l|,^,^, under some specific assump- 
tions. The assumption atlg = was considered in 
Refs. H,^; and a^ ^ in Ref. §. In Rcf. ||, a gen- 
eralization of the Weinstein criterion for the collapse [ p^ 
was obtained, in case of three dimensions, and when there 
is only cubic term in the NLSE. 

It will be interesting to derive criteria for stability more 
generally applied, using a variational approach and by 
comparison with results obtained with numerical simu- 
lations. For this purpose, we should first obtain expres- 
sions for the mean-square-root that cxtrcmizes the total 
energy. In the variational approach, for D dimensions, 
the equations for the total energy ( p0| ) and chemical po- 
tential are given by ( |2C| ) and (|2l|). So, in case of = 3, 
the equations corresponding to (|22|), ( p3| ) and (|2j) are 

al-at-Pal-Q = 0, (27) 



... X „ 9 1 P 5a2 1 Q 



p 

fi{as) = + 2^ = 



7a| 
2 



1 

2^ 



Q_ 

2a6 



(29) 
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Here we should note that an exact expression for the 
Hamiltonian, corresponding to the VA given by Eqs. (po|), 
can be derived for the general case. This will be useful 
for the exact numerical calculations that we are going to 
perform in 3D. We consider the following scaling in the 
wave- function ip of the system p3| : 



a-^/\{C), with e=-. 

a 



Then, for the total Hamiltonian, we obtain 



where 



Aq Yo Aq z,. 



(30) 



(31) 



(32) 




FIG. 2. For A2 < (with n = 0.6, or P = -0.1692) and 
A3 = (Q = 0), as in Fig. 1, we show a variational plot of 
U{a) (upper frame), as a function of the width a; and the 
corresponding phase plane defined by (at, a). 



(33) 



By taking SH{a)/6a\a=i — 0, we obtain the characteris- 
tic equation 



Xs - (r') + -X2YS + XsZs = 



With the above equation, for the energy we have 



(34) 



H 



2X, 



(35) 



Correspondingly, for the chemical potential, we obtain 



(36) 



Interesting to observe that when the two-body interac- 
tion is zero, the exact result for the chemical poten- 
tial is twice the mean square radius, as shown above [in 
agreement with the variational result of Eq.(p9|), where 
= 2/3(r^)]. As specific cases will be considered in the 
next subsections, we will discuss this point in more detail. 
For now, we can use the specific example of Q = Oi rep- 
resented in Fig. 2, in order to illustrate the validity of a 
general criterion (also valid for Q ^ 0) that we consider. 
In this case, we have just one minimum and a maximum 
of C/(a). 



A given chirp b{t) in the wave-function corresponds 
to the square root of the kinetic energy, given by the 
Hamiltonian in Eq. (|6|). With a small initial chirp, we 
should have oscillations around the minimum of U{a). 
The square of the chirp must be larger than the difference 
between the minimum and the maximum of the effective 
potential energy C/(a) for the system to become unstable. 
So, from Eq. (|26| ) and given that i/o is conserved, we can 
derive a minimum criterion for stability. The given ki- 
netic energy (related to the chirp) must be smaller than 
the corresponding variation of the potential energy, in or- 
der to maintain stable oscillations around the minimum 
(where a = a^i). The initial chirp fop niust satisfies the 
condition 



AU{a), 



hi < 



3 

2 f/(a,2 



U{asi) 



(37) 



where as2 is the value of a corresponding to the maxi- 
mum of U{a), and Gsi corresponds to the minimum of 
U{a). The maximum initial chirp is given by 



bo,, 



2 VUias2) - Ujasi] 



asi 



(38) 



From the above equation, we should also note that the 
critical values for the number of atoms, related to the 
parameters P and Q, can be obtained by the condition 
that U{as2) U{asi). We should also observe that the 
chirp parameter is introduced in the exact wave function 
by an exponential factor, as in the VA [see (|^)]. 

The frequency of the linear collective oscillations, cj^, 
is also another relevant quantity to analyze the stability 
of the condensate. When the system becomes unstable, 



6 



it reaches the value zero. From the radial variation of the 
force around a minimum a^i, obtained from Eqs. (^) and 
(p8|), Wi, in units of the trap frequency oj, is given by 







3 lb? 


a— a si 



1 


3Q 

"si 


4 


3P 






0^1 





(39) 



One of the requirements for instability of the BEG is that 
ldl is zero, when the minimum and maximum of U dis- 
appear. So, it is convenient to analyze the critical points 
when the system becomes unstable; at such points, the 
system collapses (in case we have only one minimum) or 
one of the phases (characterized by a minimum) disap- 
pears. From Eq. (pa). 



1 fu_ 

3 da? 
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(40) 



We can only have two minima (given by a^i and 0,3) and 
one maximum (at 0,2) for negative P and positive Q. In 
this case, two critical limits are possible, corresponding 
to the situations in which the position of one of the min- 
ima, ttsi or 0,3, is equal to the position of the maximum 
(one of the phases disappears). Note that in the critical 
limit we have amin = dmax — ac- In the next two sub- 
sections, we will consider the particular cases in which 
we have only one nonlinear term present in the NLSE: 
the cubic term or the quintic term. The discussion of the 
more general cases, where both terms are present in the 
NLSE, we leave for the section IV. 



A. Case of Q = 0, with attractive two-body term 

The roots for the maximum (0,2) and the minimum 
(flsi) of U can be obtained from (|2^): 

al~a, + \P\^Q . (41) 

Combining this equation with Eq. (HQ), we have 



1 (PU_ 
3 da2" 

In the critical limit |3 



1 



(42) 



and Tic = 2V7F|Pc| = 1.8965, 



-0.53499 (43) 



where Eqs. (^l]) and ( |l3| ) were used. With the above root 
of Eq. (^) , the other roots can also be easily found nu- 
merically for any value of P. Given P, once wc obtain 



the positions of the minimum {a si) and maximum (0,2) 
of U , we can calculate the value of the maximum initial 
chirp for the collapse, using Eq. 



Hi 



2)- 



(44) 



s2 



Our results, when we have only the cubic term in the 
nonlinear effective potential {Q — 0), is shown in Figs. 
2 and 3. In Fig. 2, for P = —0.1693 (corresponding to 
n = 0.6), we represent in the lower frame the phase space 
given by at versus a. In the upper frame, we have the 
total energy, in units of N{3/2)huj, as a function of a. 




FIG. 3. Results for Q = 0, as functions of n, the reduced 
number of atoms defined in the text. In the upper frame we 
have E = nil, the total energy in units of {N/n){huj/2); in 
the middle frame, the mean-square-radius (r^), in units of 
h/{muj); and, in the lower frame, the maximum initial chirp. 
In the critical limit, ric — 1.8965 and Pc = — 0.53499, the total 
energy is given by _E = nil — ricVE — 4.2407, corresponding 
to Et = y/ZNc(fiuj/2). 

In Fig. 3, we show in the lower frame our results for 
the maximum initial chirp, 6o,m, as a function of n, for 
stable solutions, n is the reduced number of atoms, given 
by Eq. (p^. In this figure, we also represent the total en- 
ergy (upper frame) and the average of the square radius, 
(r^), as a function of n. 

The static variational expression for the maximum ini- 
tial chirp, was derived from Eq. (H) and energy conser- 
vation requirements, given by Eq. (|37|). The circles cor- 
respond to results obtained by the time-dependent VA, 
by increasing the initial conditions for the chirp till the 
system collapses. We also obtain the results for the max- 
imum initial chirp using the same Eq. (|37|), but consid- 
ering exact PDE results for the observables presented in 
it. The results, for the static case, are represented by a 
solid line; and, for the time-dependent case, by squares. 
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We have observed small differences in the absolute val- 
ues of 60, m when considering the trial bo as positive or 
negative. This fact can be qualitatively understood by a 
small dephasing of the wave that occurs near the wall of 
the oscillator, when it is bouncing back. 

For the mean-square-radius, we observe that the VA 
results approaches the exact PDE results also in the un- 
stable branch in the limit n — > 0. These results are con- 
sistent with the ones obtained in Ref. |Q, where it is 
shown that, for the trapped NLSE with cubic term, the 
unstable branch of the chemical potential diverges in the 
limit 71 — > 0, with (r^) collapsing to zero (this limit is the 
non-trapped solution of the NLSE). 



B. Case of P = with attractive three-body term 

As in the previous subsection, for a stable condensate 
the maximum initial chirp isgiven by Eq. (|3^), where 
the roots are given by Eq. (^7|). The real and positive 
roots correspond to stable solutions, such that only one 
is meaningful when Q > 0. When Q < 0, we obtain two 
real and positive roots of a, in case > Q > —1/4. These 
roots are 



Osl/s2 



1 



± 



(45) 



where the minimum of U is given by a+ (= asi) and the 
maximum by a_ (= as2)- The critical limit happens for 
af = 1/2, or = -\Qc\ = -1/4. From Eqs. (g|) and 



1 

3 da2 



± 



(46) 



and the frequency of the linear collective oscillations, ujl, 
near the minimum is given by 



I 2^1~4\Q\ 
\Q\ 



(47) 



In Fig. 4, we present the main results for this sub- 
section with P = 0. As noted in Eqs.(|2|) and @) 
(respectively, the variational and exact expressions), the 
chemical potential is just twice the mean-square-radius. 
In the plots, the VA results are shown in dashed lines 
and the exact PDE in solid lines. The variational crit- 
ical limit, Q = Qc = —1/4 (a'^ = V^), corresponds to 
n3,c = 1.2371 , where 123 is the reduced number of atoms, 
that is given by the normalization of the wave- function. 
In this case, Eq. ([l^) will give us the relations between 
IQI, IA3I and n^: 



"3 



27r 



and IQI 



(48) 



Given A3, we can obtain the corresponding critical num- 
ber of atoms: 



3^3/231/4 



(49) 



a result consistent with no limit in the number of atoms 
if A3 ^ 0. 




FIG. 4. Results for P = 0, as functions of the reduced 
number of atoms ns. The total energy is given in the upper 
frame, in units of {N/ns){'hui/2); in the middle frame we have 
the chemical potential in units of {hijj)/2; and, in the lower 
frame, the frequencies of the collective breathing mode. In 
the VA, the critical limit for stability, Q — Qc = —1/4:, corre- 
sponds to n3,c =1.2371, where E — ngUc = 2\/27i3,c = 3.499 
[corresponding to Et = 2\/2Nc{hijj/2)]. 

In Fig. 4, we also observe that dfi/dn^ < (branch 
of minima) and dji/dns > (branch of maxima) cor- 
respond, respectively, to stable and unstable solutions. 
We note that localized stable structure can only exist 
for positive /i. The results in this case, that the effec- 
tive interaction contains only the trap and quintic term, 
are in agreement with the VK criterion p2| . The ob- 
served exact relation between the chemical potential and 
the mean-square-radius (/i = 2 (?£)), in the present case 
that A2 — 0, was proved in Eq. (pq). 

In the branch with the unstable solutions, we observe 
numerically that in the limit 713 — + the wave-function 
approach the form ex exp(— r2/2)/r; with /i — s- 1 and 
(r^) ^1/2. This solution is proved analytically to be 
the irregular solution of the oscillator [Recall that, when 
A2 = and A3 = 0, in Eq.(^), the regular oscillator so- 
lution gives ^ = 3; as for the irregular solution in the 
origin we have = 1 in our energy unit {huj/2)]. So, 
for the unstable solutions, with small [Ql, the Gaussian 
ansatz (started by considering the exact regular solution 
of the oscillator in the limit 713 = 0) fails. By keeping (ar- 
tificially) the Gaussian shape for the unstable solutions, 
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the radius is forced to be zero. This explains the dis- 
crepancy between the present VA and the PDE results in 
the unstable region. The more appropriate ansatz for the 
trapped NLSE (when A2 and/or A3 is nonzero) should be 
build considering the two (regular and irregular) solutions 
of the harmonic oscillator. 

Even considering the discrepancy between VA and 
PDE results for the unstable (maxima) solutions, there 
is a reasonable agreement for the results obtained for the 
maximum initial chirp 69, m, as shown in Fig. 5. In Fig. 
6, we have the variational result for the potential energy 
as a function of a, and the corresponding phase space, as 
in Figs. 1 and 2. Static predictions on PDE are also in 
good agreement for the maximum initial chirp. In PDE 
calculations we have verified numerically small differences 
between negative and positive 60. m, as in the case ana- 
lyzed and explained in the previous subsection (related 
to Fig. 3). Using Eq. ( |3^ ) with P = 0, the maximum 
initial chirp is given by 



/ 2 f._'WM] for |Q|«1. (50) 




FIG. 6. For A3 < and A2 = (g < and P = 0), with 
ns =0.5, as in Figs. 1 and 2, we present results of the VA for 
the total potential energy [/ as a function of a, and a plot of 
the phase space of at versus a. 



IV. DYNAMICS OF BEC WITH TWO- AND 
THREE-BODY INTERACTIONS 



The above result shows that, when the three-body pa- 
rameter is negative and small, the maximum initial chirp 
is cx or (X (l/v^g), as shown in Fig. 5. 




FIG. 5. Results for the maxilnum initial chirp &o,m, con- 
sidering A3 < and A2 = 0, using the VA and exact PDE 
results. 



In this section, we analyze the cases where we have 
non-zero cubic (A2) and quintic (A3) interaction terms in 
the NLSE. 

The different possibilities, relative to the signs of the 
two and three-body interactions are studied in the next 
three sub-sections, where we exclude one of the cases in 
which the signs are the same. When Q has the same sign 
of P, the basic physical picture is not essentially altered 
in comparison to the cases already discussed, where one 
of these quantities is zero. However, it is worthwhile to 
examine the changes in the collective excitations in the 
case that A2 and A3 are positive, as such observable can 
be useful to obtain informations about possible manifes- 
tation of three-body interactions. 



A. The case of repulsive cubic and quintic 
interactions (A2 > and A3 > 0) 

From Eq. (|^, one should observe that the frequency 
of the collective oscillations will decreases, as we increase 
P, or if the position of the minimum, Qsi, decreases. As 
all quantities present in Eq. ( p9| ) are positive, this fre- 
quency in the VA has an upper limited at the value 2-\/2, 
such that a reasonable frequency of oscillations around 
the minimum will be smaller than that (see also Ref. [Q). 

The effects of the unharmonicity of oscillations can be 
taken into account by the expansion of the effective po- 
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tential near the bottom Oc on the power of deviations 
y = a — ttsi . The equation for y(t) is: 



where 



ytt = -Biy - B2y^ - B^y^ 



+ — + — + — , 

'^sl ^sl '^sl 

6 lOP 28Q 



(51) 



Br = —!r- 
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si 



,10 



(52) 



Then, applying the perturbation techniques to the the- 
ory of nonhnear oscillations [ p5| , we find the frequency of 
weakly nonlinear collective excitations: 



^NL — y B\ 



1 



3^3 
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hBl 
\2B\ 



(53) 



where tr is the amplitude of the oscillations. The correc- 
tion to linear frequency is proportional to the square of 
the amplitude of oscillations of the condensate. 

The above estimate of the nonlinear oscillations in 
BEC, when both the two and three-body terms are posi- 
tive, can be a relevant piece of information to determine 
a possible manifestation of three-body interaction in the 
condensates and its magnitude. We should observe that, 
recently, in Ref. p6t , it was reported observation of non- 
linear oscillations in BEC of a gas with rubidium atoms. 

Let us estimate the amplitude of oscillations of the 
width of the atomic cloud. The points of maximum and 
minimum width are defined by the section points of the 
line Bq = const with the potential curve U{a). For large 
width, a >> 1, we can obtain the estimate 




(54) 



where o'-i-(-) = Omaximin)- K is natural to assume that 
the asymptotic value of the width is the averaged value 



B. Case of attractive cubic and repulsive quintic 
interactions (A2 < and A3 > 0) 

In this case, there is no collapse, but the condensate can 
have up to two distinct phases when IA2I is smaller than 
a critical value. From the beginning, when A3 = 0, ?7 has 
one minimum at a = agi and a maximum at a = 0^2 • 
By considering a fixed positive (and small) A3, a second 
minimum of U {asz) appears, corresponding to a denser 
phase in the condensate. As we increase the value of |P|, 
we can reach a critical value where the minimum of the 
normal phase disappears . 

When i/o/3 < AJJ = \U{as2) - U{as)\, with = a^i 
or as3, we have small amplitude oscillations near the fixed 



point as- With ffo/3 > Ai7, the character of the oscil- 
lations changes and we have large amplitude oscillations 
, due to the motion of the effective particle between the 
wall given by the repulsive three-body interaction and 
the other curve given by the quadratic potential. With 
as close to 033 and i/o/3 < AJ72 = \U{as2) — C^(as3)|, the 
condensate oscillates with the amplitude restricted from 
below by as2- 

The position of the minimum Osi of U{a) is defined by 
the equilibrium between two and three-body interactions; 
the positions of the other extrema of U , the fixed points 
as2 and ttgs are mainly defined by the contributions of the 
two- and three-body terms. These considerations allow 
us to obtain analytical expressions for such fixed points, 
that are given by the real and positive roots of the equa- 
tion given by the first derivative of U{a). 

In order to analyze the frequencies of the linear oscil- 
lations [given by Eq. (p9|)], we return to Eq (|40|), where 
P ~ ~\P\ and Q — \Q\. The physically relevant situa- 
tions occur when |P| is a fraction of one and Q << 1. 
And of particular interest are the cases where two phases 
(corresponding to two minima) are possible, at a = a^i 
and a — 0^3. 

Let us consider, for example, the particular case ex- 
amined in Table I, that we have a one-point solution 
given by a^i = \P\ and Q — a^-^ — \P\^. This exam- 
ple is chosen not only due to convenience but because it 
is not far from a more realistic situation. Usually we scale 
the normalization of the wave-function with n that is di- 
rectly proportional to So, the three-body parameter 
A3 is related to the particular point we are considering. 
The critical limit, in this case, is given by iPcil"* = (1/8) 
(|Pci| =0.5946) and Q^i = 1/64, as shown in Table I, cor- 
responding, respectively, to = 2.1078 and = 0.021 
[see Eq. (p^j. At such critical limit the frequency of 
the linear oscillations around the minimum a^i goes to 
zero [In a practical situation, before the system reaches 
this critical limit in the normal phase, it tunnels to the 
denser phase, when both energies are equal]. With |P| 
larger than |Pc|, there is only one phase (minimum) cor- 
responding to a denser phase. As we reduce the value of 
|P|, we reach a second critical point at as2 — OsS — \Pc3\, 
where 0^3 is the position of the minimum of the denser 
phase. A given value of |P| should be in between the 
limits |Pc3| < |P| < [Pel I for the existence of two phases. 

So, for the normal phase, we note that the correspond- 
ing frequencies of the oscillations are approximately given 
by the harmonic oscillator, with w^i ^ 2 (twice the trap 
frequency) . This can be seen in the previous cases where 
one of the nonlinear terms is zero and the other is small. 

Consider the case exemplified in Fig. 1, where P « 
-0.55 (n «1.95) and Q « 0.01 (53 = 0.016), and 
the corresponding roots for the minima are ass «0.386 
and Qsi wO.726. Using Eq. (^9|), we can observe that 
the frequency of the oscillations in the denser phase, 
UJL2 = 4.493, is much larger than the frequency of the 
normal phase, ujli = 1.339. If we take more realistic pa- 
rameters as, for example, 173 ~ 0.001 and n — 1.5, using 
Eqs. ([1^) and (^7|), we obtain = 0.840 (minimum of 
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the 'gas' phase), 0^2 = 0.434 (maximum) and Oss = 0.105 
(minimum of the 'hquid' phase). The corresponding fre- 
quencies of the linear oscillations of the two phases are 
ujli = 1-73 and u!l2 = 258. So, in general we observe 
that it should be expected 



1.00 



WL2 > WLi. 



(55) 



This unequality for the frequencies was also confirmed in 
numerical PDE calculations done in Ref. [||. 

The amplitude of the nonlinear oscillations in the dense 
phase can be found from the observation that Uss << 1. 



Then the terms with 



and 



are dominant in the 



effective potential energy, and an approximate solution 
can be found: 
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2|i?o| 



(56) 



An interesting phenomenon occurs when the initial 
OjIq — dgbo is large. Then, the character of cloud os- 
cillations will change, from oscillations near a^i and 0^3 
to large amplitude oscillations defined mainly by the 
quadratic potential (see Fig. 1). Let us estimate the 
criterion for the bifurcation phenomenon in the oscilla- 
tions. Taking into account the expression for the energy 
Hq, the bifurcation point is given by 



> 



2AU 1 
3 flsi 



(57) 



The result is shown in Fig. 7, using our Gaussian VA, 
which gives us a qualitative picture of the exact results, 
that were already presented in Ref. [pi . 




FIG. 7. VA results for the total energy and for the corre- 
sponding maximum initial chirp, as functions of the reduced 
number of atoms n, in case we have attractive two-body and 
a fixed repulsive three-body interactions with 53 =0.016. 
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FIG. 8. Time evolution of the square radius with the 
full PDE in the case of phase transition for n =1.754, 
(?3 =0.016 (A2 < 0, A3 > 0), beginning in the gas (less dense) 
phase solution. The initial conditions (chirp) considered are 
bo =0.05 (dotted), 60 =0.10 (dot-dashed), 60 =0.11 (solid 
line), bo =0.15 (long dashed). 

In Fig. 8 we show full numerical PDE calculation for 
the system beginning in the gas phase, which corresponds 
to the right minimum of the Fig. 1. Our results show that 
for a small chirp {bo < 0.10), the mean square radius os- 
cillates with small amplitude, according to its collective 
frequency. Applying a stronger chirp, near bo > 0.11, we 
observe transitions back and forth between two phases, 
characterized by different amplitudes of the oscillations. 
As bo is large enough we observe that the oscillation pat- 
terns remains with almost fixed and large amplitude. A 
similar picture happens if we begin in the denser phase, 
with different initial conditions given by the chirp. These 
results were obtained by using the exact PDE calcula- 
tion, and they are in qualitative agreement with the VA 
approach. The quantitative agreement between the VA 
and the exact PDE deviates considerably in the phase 
transition region. This should be expected, as we are far 
from the harmonic oscillator behavior, in a region where 
an improvement in the ansatz is necessary. 



C. Case of repulsive cubic and attractive quintic 
interactions (A2 > and A3 < 0) 

As both two and three-body nonlinear terms are non- 
zero and with opposite signs (with attractive three-body 
and repulsive two-body interactions), one could expect a 
behavior similar to the case that was analyzed in the pre- 
vious subsection, and also represented in Fig. 1 (where 
the two-body interaction is attractive and the three-body 
is repulsive). But contrary to such expectation, this case 
shows a different behavior for small radius, and no phase 
transition is possible. We observe that in this case the 
system can collapse, as the behavior for small radius is 
dominated by the three-body term, that is negative. This 
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is represented in Fig. 9, in a variational plot of U{a) ver- 
sus a, together with a corresponding plot for the phase 
space. 

In Fig. 10, we present the breathing mode collective 
excitations for a few values of (73 calculated in the vari- 
ational approximation. The collective excitations show 
that even for small negative a limited number of atoms 
is allowed in the condensate. Now, only a region of sta- 
tionary condensate - the denser phase - exists. Two points 
exist that extremize the energy: one, Osi, is stable; and 
the other, 0^2, is unstable. Here, the situation is similar 
to the cases that we have a single attractive nonlinear 
term (cubic or quintic); and the conditions for collapse, 
in terms of the initial chirp bo, is such that it must be 
larger than 60. m given in Eq. (p8[). 

One should note that this particular case (A2 > and 
A3 < 0) can be relevant in the analysis of experiments 
with BEG performed with atomic systems that have re- 
pulsive two-body interaction. No collapse is expected if 
a real three-body effect is not manifested, or in case the 
possible three-body effect is also repulsive. However, an 
attractive three-body effect will change this scenario, as 
the system must collapse for certain critical maximum 
number of atoms. In this perspective, the present anal- 
ysis shows that experiments with BEG may be useful in 
detecting negative three-body forces. The correspond- 
ing maximum critical number of atoms for stability of 
the condensate can also be obtained from the present ap- 
proach, by taking the limit &o,m — > (or — > flsi)- For 
instance, when 33 — —0.5, the critical number n can be 
obtained from the rhs (positive A2) of Fig. 11. 
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FIG. 10. Collective excitations for the breathing mode of 
the condensate for positive two-body interaction (A2 > 0), 
and for a set of negative three-body interaction (A3 < 0), as 
shown inside the figure. 
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FIG. 11. Total energy E, mean-square-radius (r^), and 
the chemical potential /.i, for negative A3 (with 33 = —0.5). 
The dashed line corresponds to the variational approach and 
the solid line to exact numerical calculations. Departing from 
/I =3, to the right or to the left, we have the solutions corre- 
sponding to minima of the total energy, until critical values 
Tic (in both sides) are reached. After the critical numbers, the 
curves follow lower branches, corresponding to maxima. 



FIG. 9. As in Figs. 1, 2 and 5, we have in the upper 
frame the potential energy (7(a); and, in the lower frame, a 
plot of the corresponding phase space of at versus a. We have 
nonzero values for both two and three-body parameters, with 
A2 > and A3 < (53 = -0.5). 



We present in Fig. 11, for a fixed attractive three- 
body interaction, with (73 — —0.5, the total energy _E, the 
mean-square-radius (r^), and the chemical potential, as a 
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function of the reduced number of atoms n multiplied by 
the sign of the two-body interaction. We consider both 
negative (A2 < 0) and positive (A2 > 0) two-body in- 
teractions, and the value of 173 = —0.5 was chosen due to 
numerical convenience. The variational approach (dashed 
line) gives a good description of the minima for small |n|, 
but fails particularly for the maxima solutions, as one can 
see in the figure when comparing with the exact results 
(solid lines). In agreement with the variational results 
of collective excitations, the number of atoms allowed is 
limited to a critical number 3.5 (the corresponding 
exact result is ~2.93). For a more realistic value of the 
three-body parameter, with 173 = —0.01, the exact critical 
number is ric ^^4200. 

In the upper frame of the Fig. 11, it is shown a sim- 
ilar pattern for the minimum and the maximum of the 
energy, as already described with A3 > when A2 < 
and with A3 < when A2 = 0. Such results also indicates 
that an initial wave function with enough larger chirp can 
make the condensate unstable and that the correspond- 
ing magnitude of 60, m can be inferred from the absolute 
diff'erence between the maximum and minimum of the en- 
ergy. This implies that, when we have positive two-body 
interaction (A2 > 0), the presence of a negative three- 
body term (A3 < 0) can be detected in principle even if 
the number of atoms has not achieved the critical num- 
ber. The initial maximum chirp plays a relevant role in 
this case. 



V. VK CRITERION 

In the context of stability analysis of NLSE bound 
states, the orbital stability, or stability with respect to 
the form, was first settled by the so-called Vakhitov- 
Kolokolov criterion that is detailed in a recent re- 
view in Ref. |l^. A solitary wave is orbitally stable if 
the initial orbit, chosen near the ground-state orbit, im- 
plies that the orbit of the solution at any t > remains 
close to the ground state. The criterion was first demon- 
strated in Ref. |33|, by minimizing the deviations between 
the orbital states and the ground-state, with respect to 
the initial position parameters. In our notation, the VK 
criterion for stability can be written as ^ 



dN 



< 0. 



(58) 



However, the original VK criterion and the above cited 
demonstration have only considered cases without trap, 
with nonlinearities expressed by /(IV'P)- More recently, 
the general conditions that guarantees orbital stability of 
stationary trapped condensates, described by the NLSE 
with cubic term, was formally demonstrated by Berge et 
al. ||2^], in agreement with results analyzed in Ref. j^] for 
the case that the cubic term is negative. 



The validity of the VK criterion appears to be solidly 
consolidated by the above described formal demonstra- 
tions and by the specific cases that have so far been con- 
sidered. These results appear to support an extended 
range of applicability of the VK criterion to the NLSE. 
In the several cases that we have considered in the present 
paper (by exact numerical procedure and also by the 
Gaussian variational approach), the applicability of the 
VK criterion is out of the question in the mentioned cases 
that have been examined by other authors. We have also 
confirmed the applicability of the criterion in the case 
that we have only the trap and the quintic term in the 
effective interaction, as shown in Fig. 4. 

The validity of the VK criterion cannot be extended 
to trapped NLSE with two and three-body terms, when 
the two body term is negative, as one can observe from 
the results obtained in Ref. . This is an unexpected re- 
sult, considering that the criterion is applicable for non- 
trapped systems (with cubic and/or quintic terms) and 
also for the trapped cases when we have only the cubic or 
the quintic terms (and also in the cases that both terms 
have the same sign) . The VK criterion is also not applica- 
ble when A2 > 0, as observed in Fig. 11, in the right side 
of the plot. The criterion fails particularly in the region 
near the oscillator solution (where fj, =3), where we can 
observe that dfi/dn > in a stable branch (upper) and 
also in the unstable branch (lower). (We recognize the 
stable branch as the one that corresponds to minima of 
the energy.) Therefore, in the present work, we confirm 
the conclusion about the limitation of the VK criterion: 
The nonvalidity of the VK criterion is verified numeri- 
cally when we have a positive harmonic trap, with two 
and three-body nonlinear terms with opposite signs. 



VI. CONCLUSIONS 



The stability of a trapped condensate, with two- and 
three-body nonlinear terms, was studied in the present 
work, considering several aspects, as the initial conditions 
in the wave-function and the validity of the Vakhitov- 
Kolokolov criterion. For the initial conditions we consid- 
ered a chirp parameter, which is related to the initial fo- 
cusing (defocusing) of the cloud. A non-zero initial chirp 
bo introduces oscillations in the condensate near the min- 
imum of the energy, such that it can lead a previously 
stable system to collapse, when bo reaches a maximum 
limit 60, m, that corresponds to the energy difference be- 
tween a minimum and a maximum of the total potential 
energy. So, in the presence of a fixed value of the chirp, 
the number of atoms of a stable condensate with attrac- 
tive interaction is smaller than the corresponding value 
when the chirp is zero. In case that the potential has 
two minima (when cubic and quintic terms are present 



^We should note that fi corresponds to — A of Ref. 112 
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in the NLSE), related to two phases of the condensate, a 
bifurcation phenomenon is predicted in the dependence 
of the value of the initial chirp. A maximum initial chirp 
can affect the system in such a way that the oscillations 
can switch from the gas phase to the liquid phase and 
vice- versa. 

The present study was performed by using exact nu- 
merical solutions of partial differential equation, as well 
as by a corresponding variational approach. The analyt- 
ical predictions based on the time-dependent VA, using 
Gaussian ansatz, are qualitatively confirmed by the exact 
time-dependent numerical simulations. We would like to 
point out the advantages of the Gaussian VA when com- 
paring with the moments method, which encounters some 
difficulties in deriving the collapse conditions even in the 
pure attractive two-body case, for some classes of the ini- 
tial data (the sign of the time dependent function in the 
equation for the mean-square-radius, {r"^), is not defined 



The VA starts to deviate from the exact results when 
the system is close to the collapse conditions and partic- 
ularly for the unstable solutions (maxima) of the total 
energy. For this regions, where the Gaussian VA fails, 
one should improve the ansatz or take the results as a 
qualitative picture to guide the exact numerical calcu- 
lations. In this respect, one should note that, for the 
NLSE without the trap potential, it is well known that 
the time-dependent variational approach fails to describe 
the region near the collapse . 

Considering that in this work we have shown the rele- 
vance of an initial chirp parameter to study the stability 
of the condensate, we should observe that a good esti- 
mate for such parameter relies in a good approach to 
determining both minimum (stable) and maximum (un- 
stable) solutions for the total energy. The Gaussian VA 
has a good agreement (quantitatively) with exact results 
for the minima, but can only give a qualitative description 
for the maxima. Still we can observe that the VA calcula- 
tions of the chirp parameter, when compared with exact 
numerical calculations, are in reasonable agreement. 

We also would like to point out that the present ap- 
proach can be extended to study the stability of chirped 
laser beams in inhomogeneous three dimensional media 
with Kerr nonlinearity. 

Finally, the main results of the present work are: (i) A 
chirp parameter in the wave-function is shown to be use- 
ful to study the initial conditions for a stable condensate 
to remain stable, (ii) we show that the VK criterion can- 
not be extended to cases of harmonic trapped BEG when 
the nonlinear two and three-body terms have opposite 
signs; (iii) if an atomic system has repulsive two-body 
interaction, the collapse is possible if the effective three- 
body interaction is negative. In this perspective, one can 
use the observed critical number of atoms in order to de- 
termine the corresponding three-body parameter. 
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